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Abstract 

Perturbation of gravitational fields may be decomposed into scalar, vector and tensor compo- 
nents. In this paper we concern with the evolution of tensor mode perturbations in a spatially 
closed de-Sitter background of Robertson- Walker form. It may be thought as gravitational waves 



in a classical description. The chosen background has the advantage of to be maximally extended 
and symmetric. The spatially flat models commonly emerge from inflationary scenarios and are 
not completely extended. We first derive the general weak field equations. Then the form of the 
field equations in two special cases, planar and spherical waves, are obtained and their solutions 
are presented. We conclude with discussing the significance of the results and their implications. 
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I. INTRODUCTION 



Here we investigate the freely propagating gravitational field requiring no local sources for 
their existence in a particular background. As an essential feature of the analysis of general 
theory of small fluctuations, we assume that all departures from homogeneity and isotropy 
are small, so that they can be treated as first order perturbations. We focus our analysis 
on an unperturbed metric that has maximal extension and symmetry by taking K = 1 and 
presence of a positive cosmo logical constant. The background is de-Sitter spacetime in slicing 
such that the spatial section is a 3-sphere. In the previous works mostly the case K — wre 
investigated extensively [1-3]. Even though in some works K is not fixed for demonstrating 
the general field equations, but for solving them usually K = is imposed [4,5]. The 
study of this particular problem is interesting and relevant to present day cosmology for the 
following. Seven-year data from WMAP with imposed astrophysical data put constraints on 
the basic parameters of cosmological models. The dark energy equation of state parameter is 
— 1. 1 ±0. 14, consistent with the cosmological constant value of —1. While WMAP data alone 
can not constraint the spatial curvature parameter of the observable universe Qk very well, 
combining the WMAP data with other distance indicators such as H , BAO, or supernovae 
can constraint Qk- Assuming u = —1, we find il\ = 0.73+/— 0.04 and Vt to tai = 1.02+/— 0.02. 
Even though in WMAP seven-year data it has been concluded as an evidence in the support 
of flat universe, but in no way the data does not role out the case of K = 1 [6]. Planck 
data report release in near future makes the situation more promising to settle down this 
dispute. In the analysis of gravitational waves commonly Minkowski metric is taken as 
the unperturbed background. According to mentioned observational data, the universe is 
cosmological constant dominated at our era. So in the analysis of gravitational waves we 
should replace the Minkowski background with de-Sitter metric. The essential point is 
that spatially open and flat de-Sitter spacetime are subspaces of spatially closed de-Sitter 
space. The first two are geodesically incomplete while the third is geodesically complete and 
maximally extended. From the singularity point of view the issue of completeness is crucial 
for a spacetime to be non-singular. Taking the issue of completeness seriously, we have no 
way except to choose K — 1. By choosing the maximally extended de-Sitter metric as our 
unperturbed background we include both cosmological and curvature terms in discussion of 
gravitaional waves [7-11]. We begin by deriving the required linear field equations. 
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II. LINEAR WEAK FIELD EQUATIONS 



Supposed unperturbed metric components in Cartesian coordinate system are 12] : 



9oo = -1, 9io = 0, gn = a {t)g. 



Kl 



'X. 'X'.' 



a(t) = acosh(t/a), g {j = + if - - 2 , (1) 

with the inverse metric 

0°° = -l, ^ 0i = 0, ^' = «W, g? = {&-K!te) i (2) 

where if is curvature constant and a = The non-zero components of the metric 

compatible connections are: 



r°. = aa(8ij + K — ) = ad&j 



p - -5- 

rj fc = fj fc = (3) 

Dot stands for derivative with respect to time. Since we are working in a holonomic ba- 
sis, then the connection is torsion-free or symmetric with respect to lower indices. Let us 
decompose the perturbed metric as: 

9nv = 9vv + h pu , (4) 

where g^ is defined by Eq.(l) and h pu is small symmetric perturbation term. The inverse 
metric is perturbed by 

x ^ = gr - gr = -r p g va h pa , (5) 

with components 

X 00 = -hoo, x m = a 2 (hm-Kx i x i h j0 ) 

x v = -a-\hij - Kx i x k h k] - Kx j x k h ki + K 2 x i x j x k x l h H ). (6) 



Perturbation of the metric produces a perturbation to the affme connection JjJ 



$Kx = ^"1-2VCa + <9aV + - d p h Xu \. (7) 
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Thus Eq.(7) gives the components of the perturbed affine connection as 

"| ^3 3 ^ 

5T) k = — [-2aa(h i0 - Kx l x l h l0 )(5 jk + K JZ^ ) ~ 2K ( h *m - Kx l x l h lm )(5 jk + K - % _ \ x2 )x m 



+d k hij + djh ik - dihjk - Kx l x J (d k hij + djh ik - dih jk )} (8) 
Sr j0 = ^(~ 2 ^ h ij + hij + djh i0 - dih j0 + 2Kx i x k -h ki - Kx i x k h kj - Kx i x k d j h k0 + Kx i x k d k Ug) 

SFl = -{2aa(5 tJ + X — j^L_ )fr 00 + 2Kh ok (S ij x k + K *** x2 ) ~ d t h 0l - djh 0i + h l3 \ (10) 

5r 00 = ^ [2h 0i - dihoo - 2Kx i x j h 0j + KxWdjhoo] (11) 

5T% = -h i0 - Idihoo (12) 

Ko = —hoo (13) 
The tensor mode perturbation to the metric can be put in the form 

h 00 = 0, h i0 = 0, hij = a 2 D i:j , (14) 
where D^s are functions of X and t, satisfying the conditions 

g'Wij = 0, <,' l V,I) jk 0. (15) 
The perturbation to the affine connection in tensor mode are 

5T° 00 = 5T% = 5V 00 = (16) 
Sr% = aaD l3 + yAj (17) 

5T) = aaD i:j + yAj (18) 
<Sr$ fc = \\dkDij + djD ik - diD jk - 2K(D im - Kx l x l D lm ) x 

(5 jk + K x X ^ x2 )x m - Kx i x\d k Dij + djD lk - d t D jk )]. (19) 
The Einstein field equation without matter source for the tensor mode of perturbation gives 

5R jk = -Aa 2 D jk , (20) 



where 



3 • a 2 ■■ 
SR jk = -(2d 2 + ad)D jk - -D jk - —D jk 

+ l -O l d l D jk - 4KD jk - y {didmD^xW- 

~Kx m d m D jk - K(d k D mj + d 3 D mk )x m 

x „mJL 



+K 2 D ml (S jk + K- ^)x m x l . (21) 

1 — Kx z 



Scale factor a(t) satisfies the Friedmann equation, so we get 



2d 2 + ad = Aa 2 - 2K. (22) 
Inserting Eq.(22) in Eq.(21) and Eq.(21) in Eq.(20) , we would have 

~aaD jk - yi), fe - 2KD jk + ^&diD jk - ^{d i d m D jk )x i x m - %Kx m d m D, jk 

j k 

-K(d k D mj + d j D mk )x m + KD ml (5 jk + K l X _ X Kx2 > m ^ 1 = (23) 

It is straightforward to show that 

^V 2 D jk = ^V m V n D jk = ^d t dW ]k -j(d l d m D jk )x i x m -KD jk 

-Kx\d k D l3 + djD lk ) - ^Kx'diDj, + Kx l x\5 jk + K ^ 2 )Ai- (24) 

It remains to put Eq.(24) in Eq.(23), then we obtain the final equation. 

1 3 a 2 ■■ 

-V 2 D jk - -aaD jk - —D jk - KD jk = 0. (25) 

Our first task to establish the field equations is fulfilled. Next we look for special solutions 
of this field equation analogue to plane and spherical waves. For the plane wave like so- 
lutions, using Cartesian coordinate systems is suitable while for the spherical waves, polar 
coordinates (%, 9, 0) are convenient. 



III. PLANE- WAVE ANALOGUE 



In the case of flat models i.e. K = condition (15) reduces to 



Ai = 0, 6 ik diD kj = 0. 



(26) 



Looking for a wave propagating in z-direction, Eq.(26) simply gives 



A 3 = o, 



(27) 



with two independent modes 



D ij+ = D + (z,t) 



^10 0^ 

0-10 
y0 Oy 



, and D ijx = D x (z,t) 



^010^ 

1 
yO Oy 



(28) 
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where D(z,t) satisfies D 2 D(z,t) = with the well-known plane wave solution. 

In the case of K — 1 an analogue solution for Eq.(15) exists. Following a lengthy calculation 

due to non-diagonal components of ^ we obtain 

^ l n xz ^ 

i — zn — it U 



Vi -x 2 

where X 2 = x 2 + y 2 + z 2 and 

D x (z,t) 



l—x' 2 —z 2 





(l-z 2 )(l-x 2 -z 2 ) 



l—y 2 —z 2 



-yz 



(l-Z 2 )(l-y 2 ~Z 2 ) 

_ yz z 2 (x 2 -y 2 ) 

\ (l-z 2 )(l-y 2 -z 2 ) (l-z 2 )(l-y 2 -z 2 ) (l~z 2 )(l~x 2 -z 2 )(l-y 2 ~z 2 ) ) 



(29) 



D 



XI] 



VI - X 2 (l - y 2 - z 2 ) 




1 



1 

2xy 



yz 
T^T 2 

xz(l+y 2 —z 2 ) 



\ 



l-y' 2 -z' 2 (l-z' 2 )(l-y' 2 -z 2 ) 
yz xz(l+y — z ) Ixyz 2 

\ T^z 2 (l-z 2 )(l-y 2 -z 2 ) (l-z 2 )(l~y 2 -z 2 ) / 



(30) 



By inserting Eqs.(29) and (30) in Eq.(25) with some manipulation we conclude that each 
mode, x and +, satisfies 



( 1 - Z )Q^ D ^ t )+^Q- z D (^ t )- D ^ t )+ | _ .2 

-3ad J D(z,t) - a 2 D{z,t) - 2D(z,t) = 0. 



(31) 



We use method of separation of variables to find the solutions of Eq.(31). Then we may 
write 

D(z,t) = D(z)D(t) (32) 

Using Eq.(32), Eq.(31) leads to 

1 - z 2 d 2 



D(z) dz 2 



D(z) + -^--?-D(z)-l + 



D(z) dz 



1-z 2 



2 D(t) D(t) 
a 2 -M + 3aa^M + 2 



D(t) D(t) 

Eq.(33) may hold merely if each side is equal to a constant, i.e. we have: 

D„(z) + — - 1 + = -q 2 



D q {z) dz 



D q {z)dz 2 q 



a 2 D(t) 3aa ^ . . 9 
+ ^— A,(t) + 2 = -g 2 , 



1 - 



(33) 



(34) 



(35) 



D q (t) D q (t) 

where q 2 is an arbitrary positive constant. We should take it positive since we are looking 
for a periodic wave. Eqs.(34) and (35) can be written as 

d 2 „ , x „ 9 „ , x , „ 6 



(36) 



and 

a 2 D q (t) + 3aaD(t) + (q 2 + 2)D(t) = (37) 
To solve Eq.(36) and finding D q (z) we recognize that it can be given by 

D q (z) = (l-z 2 )-^U q (z), (38) 

where U q (z) are solutions of the following equation 

(1 - z 2 )^U q (z) + *z±U q (z) + (q 2 - A)U q {z) = (39) 

Let us define V q (z) = -^U q {z) and differentiate Eq.(39) two times, this gives 

(1 - ^)^V q (z) - zj-V q {z) + q 2 V q (z) = (40) 

It should be noticed that Eq.(40) satisfies Eq.(36) provided that we replace q 2 by q 2 — 3 in 
Eq.(36). We recognize that the solution of Eq.(40) for V q (z) can be written as Chebyshev 
polynomial, provided that q accepts integers, i.e. q = n, so that 

V n (z) oc exp(±m arccos z). (41) 

Then we have 

D n (z) = (1 — z 2 ) exp(±in axccos z) . (42) 

Next we examine the temporal function of the mode. Imposing the condition that q 2 replaced 
with q 2 — 3 and q be integer, gives 

a 2 D n (t) + 3aaD n {t) + {n 2 - l)D n (t) = 0. (43) 

It is convenient to define the conformal time r by: 

dr = — ^-y where a(t) = acosh(t/a;) (44) 

Integrating Eq.(44) gives 

exp(t/a) = tan(r/2) (45) 

Notice that t = — oo, 0, +oo corresponds to r = 0, n/2, n respectively. So while the domain 
of time is — oo < t < +oo, the domain of conformal time is < r < ir. 
We may recast Eq.( 43) in terms of conformal time as 

^D n (r) - 2 cot r^D n (r) + (n 2 - l)D n (r) = 0, (46) 



where D n (t) = D n (r). To solve Eq.(46), let us define a new parameter Y = cosr with 
domain — 1 < Y < +1 and t = — oo, 0, +oo correspond to Y — 1,0,-1 respectively. In 
terms of the new parameter Y, Eq.(46) becomes 

(1 " Y2) w^ 3n{Y) + Y w Dn{Y) + (n2 " l)bn{Y) = °' (47) 

where D n (r) = D n (Y). 

If we define W n (Y) = -^D n (Y) and differentiate Eq.(47) with respect to Y, this gives 

(! - y2 )^~2 W ^ Y ) ~ Y ^ W ^ + ^(Y) = (48) 
Again Eq.(48) is Chebyshef of first kind and its solutions are 

W n (Y) = exp(±in arccos Y) (49) 

For last step we should solve 

-^—D n (Y) — exp(±in arccos Y) (50) 

(Jj 1 

We have previously defined Y = cosr, so we have r = arccos Y. Let us recast Eq.(50) in 
terms of r, it becomes 

±D n (r) = -smre ± ^, (51) 



with 



We may write 



DJt) = — ^t(cost =f insmr)e ±inT . (52) 
1 — n 2 



~ (l-z 2 ), . I exp[±in(aiccos z + t)] 

D n (z,T) = — — (cosr ± insmr) < (53) 



1 -n 2 



I exp[± 



in arccos z — r 



The first mode n — 1 is pure gauge mode and should be excluded from the acceptable 
solutions. This is the analogue of a plane wave moving in z-direction for a closed model. 
We may find a similar solution for the waves that are analogue to plane wave in x and y 
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directions. In this case we would have 



/ 



D+(x,y,z,t) 



Vi -x 2 



x 2 (z 2 ~y 2 ) — xy xz 

(l-x 2 )(l-x 2 -y 2 )(l-x 2 -z 2 ) (l-x 2 )(l-x 2 -z 2 ) (l-x 2 )(l-x 2 -z 2 ) 
-xy -1 



(l-x 2 )(l-x 2 -z 2 ) 



\ (l-x 2 )(l-: 
( 



(l-x 2 -y 2 ) 









1 



(54) 



D x (x,y,z,t) 



D x {x,t) 



(1 -x 2 -y 2 )Vl -X 2 
( 



2x 2 yz 



(l-x 2 -z 2 ) J 
xz(l-x 2 +y 2 ) xy \ 



(l-x 2 )(l-x 2 -y 2 ) (l-x 2 ){l-x 2 -y 2 ) 1-x 2 



xz{l—x z +y 2 ) 
(l-x 2 )(l-x 2 -y 2 ) (l-x 2 -y 2 ) 



2yz 



xy 

\ TT^x 1 ) 



D + (x,y,z,t) = 



D+(y,t) 
Vi -x 2 



1 

1 — x 2 — y 2 
xy 



{1 - y 2 ){1 _ x 2_y2 ) 

y 2 (x 2 -z 2 ) 



1 




-yz 



(55) 



V 



(l- 2/ 2 )(1 _ 2/ 2_ 2 2) (l-y2 ){l ^ y 2_ z 2 ){l _ x 2_y2 ) {1 _ y 2 ){1 _ y 2_ z 2 ) 

o - 1 



(l-y 2 )(l-y 2 -z 2 ) 

( - 



,(56) 



1 — y — z 



D x (x,y,z,t) = 



D x (y,t) 



(1 - y 2 - z 2 )Vl - X 2 
where D(x,t) and D(y,t) are given by: 



l-y 2 



1 



\ 



yz 2xy 2 z xy(\-y 2 -z 2 ) 

T^F (l-y 2 )(l-y 2 -z 2 ) 

i xy(l-y 2 +z 2 ) 2xz 

\ 1 (l- y 2 )(1 _ 2/ 2_ z 2) l-j,2_ z 2 y 



(57) 



and 



1 - x 2 

D n (x, t) = -(cost ± in sin r) 

1 TX 



l-y 2 

D n (y,t) = - r(cosr ± in sin t) 

1 — n z 



{exp[±in(arccosx + r)] 
exp[±in(arccos:r — r)] 

{exj9[±in(arccosy + r)] 
exp[±m(arccos y — r)] 



n ^ 1 



(58) 



(59) 



It is important to notice that Eq.(54) and Eq.(55) may be achieved from Eq.(29) and Eq.(30) 
respectively by the coordinate transformation x — > z, y — >■ y, and z — > —x. This leads us 
to write the solution of waves moving in an arbitrary direction. Let us assume this arbitrary 
direction is 

hi = sin 9 sinip, fi2 = sin 9 cosip, fi3 = cos 9 (60) 

The result can be obtained from Eqs.(29) and (30) by the coordinate transformation intro- 
duced by 

z ->■ n • X 

hz{h ■ X — z 



y 



X ->■ 



5,2 



n 2 a; - niy 



(61) 
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We get 

fl — (h-X) 2 ) exp[±m(arccos(n • X) + t\ 

D n+ (X,t)=A+ {X,n) y { ^' ' (cos r±m sin t)x\ n V V J n ^ 1, 

Xi i Xi i In I exp[±m(arccos(n • X) - r] 

(62) 

— * 

where the explicit forms of A+ (X, n) are listed in the appendix. 

Xjj 

This result may be used to expand a general function as linear superposition of these eigen- 
functions, i.e. it should be replaced with exp(ik^x^) in Fourier transformations. 



IV. SPHERICAL WAVE ANALOGUE 

To consider this case it is more suitable to work in polar coordinates, x 1 = (%, 9,(f>). In 
this basis the non-zero components of the unperturbed metric are: 



9u = 1, 922 = sin z x, 933 = sin 2 x sin 2 9, 



(63) 



with the inverse 

-11 1 -22 • 2 -33 • -2 ■ -2 a 

9 = 1, 5 =sm X, £ =sm xsm 0. 
The non-zero components of the unperturbed connections are 

T\ 2 = - sin x cos x, r^ 3 = - sin x cos x sin 2 0, 
r|i — cot x, T| 3 = — sin 9 cos 0, 

r! 1 = cotx, ri 2 = cot0. 

In this case g^ is diagonal and the conditions (15) for a transverse wave give 

D u = 



(64) 



(65) 



(66) 



We may distinguish two independent polarizations as 



^0 ^ 



D +ij ( x ,9,t) = 



D xij ( X ,9,t) 



sin 2 9 



sin 9 



1 
-sin 2 



/(D 0^ 
1 

V 010 / 



(67) 



(68) 
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Inserting Eqs.(64) and (65) in Eq.(15) and expressing V 2 in polar coordinates, with a rather 
lengthy but straightforward calculation it can be shown that both D+(x,t) and D x (x,t) 
must satisfy the same equation as 

^D( X , t) - 2 cot x^-D( X , t) + 2^%^ - SaaDfa ~ a?D( X , = (69) 
d X 2 ox sm X 

To solve Eq.(66) we may assume that 

D( X ,t) = D( X )D{t) (70) 

Then we have 

1 d 2 cot X d 2 2 D(t) .b{t) 

-Cix) - 2-r--^- J D(x) + — ^— = a 7— + 3aa-^— . (71) 



D( X )dx 2 D( X )dx sin 2 x D(t) D(t) 

Eq.(68) holds provided that each side is equal to a constant, i.e 

1 d 2 O cotx d 2 



^(X) " 2^7tW^(x) + — ^- = "Of - 1)> (72) 



D(x)9x 2 D(x)dx sin 2 x 

q^( f ) 3 a ( f )a W ^ 

£>(«) £>(<) ' 

where q 2 is an arbitrary positive constant. So we have 

d 2 „, , _ 9 „, , , „ 2 



2J D( X )-2cot X ^(x) + (r + -^)^(x)=0 (74) 
a 2 (t)l)(t) + 3a(t)a(t)D + (q 2 - l)D(t) = (75) 

To solve Eq.(71) for D q (x) we may define a new parameter X = cos X and D(x) — D(X), 
then we have 

(1 " X 2 )^- 2 D(X) + 1^(1) + (g 2 - 1 + Y^)D{X) = (76) 
Eq.(71) has a solution 

D q {X) = {l-X 2 f/ 2 ±U q {X) ) (77) 

where U q (X) are solutions of 

(1 " X 2 )-^U q (X) + X^U q {X) + (g 2 - l)U q (X) = (78) 
Now V q (X) = £U q (X) satisfies 

(1 " X 2 )^V q {X) - X^V q (X) + (g 2 - l)K(X) = (79) 
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Eq.(76) is a Chebyshef type I provided that we take q = n. Then we have 



V n (x) — exp(±in arccos X) (80) 

and 

D n (x) = smxexp(±inx) (81) 
The temporal part is the same as plane wave analogue Eq.(52) and we have 

(cos r± in sin r) . J exp{±in{ X + t)) 
D n {x,t) = sinx (82) 



1 -n 



exp(±in(x - r)) 



It is interesting to note that in the case of flat models, i.e. K — 0, Eq.(71) in the (r,9,4>) 
bais takes the form 

^D(r)-l±D(r) + lD(r) = -<?D(r), (83) 

which has the solution 

D q (r) oc re ±iqr , (84) 
where q can be any arbitrary real number. If we consider the ratio ^ we get 

^ J_{-K±-(x + r)) far If = 1 (85) 
522 smx exp(±m(x-r)) 

^22 1 

— oc -exp(±iqr), for if = 0. (86) 
922 r 

Both Eqs.(82) and (83) are converged for large distances. But these solutions are not 
physically acceptable, because they are forbidden according to general theory of relativity. 



V. DISCUSSION 



Our investigations show that in analysis of gravitational waves the background of de- 
Sitter with K = +1 fundamentally differs from the scale- free de-Sitter with K = 0. We 
found the wave numbers should be discrete as already has been realized that the spectrum 
of the Laplacian in spherical space is always discrete [13]. Another relevant feature is the 
existence of cut off on the long-wavelength of gravitational waves. This may be tested by 
the measurement of dipole and higher multi-pole moments of the CMBR anisotropy which 
contains information about the long-wavelength portion of the spectrum of energy density 
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produced the large scale galactic structure of the universe. These are sensitive to the presence 
of long-wavelength perturbations. The obtained eigenmodes are the fundamental tools of 
analysis of cosmic evolution of perturbations in spatially closed models. In the formation 
of large scale structure and study of anisotropics of CMBR we should use these eigenmodes 
to expand the perturbations. Another feature of the obtained eigenmodes is that they 
are effectively transverse in zone near to origin and at far distances, in contrast to flat 
spacetimes, this is not so. For example a wave moving in z direction its 3j components are 
not vanishing exactly but in near zone may be approximated to zero. It has been shown that 
hij = a 2 Dij oc D "^zf ■ This means that in a collapsing phase i.e. in the time interval — oo < 
t < corresponding to conformal time < r < n/2, the perturbation is decaying while in 
the expanding phase, i.e. in the time interval < t < +oo corresponding to conformal time 
7r/2<T<7r,itis growing. Of course always the smallness conditions of the perturbation 
with respect to the unperturbed metric holds. The amplitude of the gravitational wave 

. / l+(n 2 — 1) sin 2 r 

changes with time as ha oc v , , ^ . 2 . So its growth is significant for the modes that 

o hj \n — 1} sin t ° ° 

satisfy the condition \/n 2 — 1 sin 2 r 1, where it changes as \hij\ oc ( - w2 _ 1 ) sin 2 T - For those 
modes that satisfy the condition y/n 2 — 1 sin 2 r ^> l,we have \ha\ oc ■ — and changes 
are relatively smooth. Singularities appear in the solutions are of coordinate type, where 
we may see them in the unperturbed metric too. At last the obtained results are crucial to 
expand any perturbation appears in different contexts of the spatially closed cosmological 
models in terms of them. 
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Appendix 



The components of the amplitude of gravitational waves moving in an arbitrary direction: 



"2 * 2 ~ 2 

— * — * H —* Tl Tl 

A + (X,h) = D + (X,n)-^ 1 + D + (A,n)-^ 

xil xii 1 — ni x 2 2 1 — tl 



D + (X,n)nl + 2D + (A,n)^?f 
X33 X12 1 - ni 



u 3 

2D + (X,h)^t= + 2D + (X,n)^= (A.l) 



X13 ' y/l -fl 2 3 X23 ' y/l-fl% 



A + (X,h) = D + (X,h)-^ + D + (X,h)^- 2 

X22 xil 1 — ^3 x 22 1 — Ti\ 

+D + (X,h)hl-2D + 

X33 X12 I—Tig 

-2B+ (X,ft)^2^ + 2D + (.?,n)^Mi= (A.2) 

X13 V 1_Tl 3 23 V 1 -^! 

A + (X,n) = D + (A,n)(l-n 2 ) + J D+ (X,n)n 2 

x 33 x 22 x 33 



-2D+ (X,n)n 3 Jl -hi (A.3) 

X23 V 

„ A. ~ ~ A. 2 

A + (X,n) = -£>+ (X,n)-^ f + J D + (X,n)^5 + J D + (X,n)nira 2 

X12 xil 1 — Til X22 1 — Tl\ X33 

+D+ (X,h)(^$)h 3 + D + (X,h)(^f£ 2 
xi2 V 1 -^/ X1 3 ^1 -n§ 

+2 fl+ (X.ft)^ ' (A.4) 

X23 s/l-nl 

A+ (X,n) = —D+ (X,h)hih 3 + D+ (X,h)hih 3 — D+ (A,n)n 2 

X13 X22 X33 X12 

D. (X,nA + D + (i,ft)*lMzl) (A.5) 

X13 x 23 V 1_Tl 3 

(A,n) = -_D+ (A, h)h 2 h 3 + £>+ (X,h)h 2 h 3 + £>+ (A,n)ni 

X23 X22 X33 x 12 

- D+ (.Y,n)4i, + D+ (X,n)^fci> (A.6) 

xi3 v 1 -™! X23 v 1 -™! 
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where -Dh 



D +11 (X,n 
D xl2 (X,n 
D +13 (X,h 

D xl3 (X,n 

D +22 (X,n 
D X22 (X,n 

D+ 23 {X,n 

D X23 (X,n 
D + 33 {X,h 

D X33 (X,n) 



X, h) are as follows: 



Vi -x 2 (i 



{h 2 x - hxyf 
l-nl 



[h-Xf) 



{VT^ll - ^■X}-^ - { n ■ X) 2 ]}- 1 



l-n§ 



(n ■ X)(n 2 x - hiy) 



y/T=X*y/T=mi - (n ■ - - (n ■ X) 2 ] 

{h ■ X){n 3 {h ■ X) - z) 

VT^VT^mi - iM t^[ Z)2 - (n • - (n ■ Xf) 

-{VT^[1 ~ ^;X}- Z)2 - (n ■ xf]}- 1 

1 n 3 

2(n 2 x - h 1 y)(n 3 (n ■ X) - z) 
(1 - hl)y/T=X*[\ - ( " 3( tl^ )2 - (n ■ X) 2 } 2 
(n 3 (n-X)-z)(n-X) 



, D xll (X,n)=0, D +12 {X,n\M§ 

(A.8) 
(A.9) 

(A.10) 



(h 2 x-h 1 y) 



n-X)[l + 



_ (n 3 {h-X)-z)* 
(h 3 (h-X)-z) 2 V,-; 



VT^(1 - fi§)(l - (n • X) 2 )[l - ^f^p - (n • X)2] 



- (n-X) 1 



vTTx*(i - (fi • x) 2 )[i - _ (a . 

(n ■ X) 2 [{n 2 x - h iy f - (n 3 (n • X) - z) 2 }{{\ - n 2 )Vl-X 2 (l - (n • X) 2 ) 

x [i - {h :r h \f - (- • - - (« • x) 2 ]}- 1 

(1 — n 3 ) 1 — n 3 

2(n 2 x - hiy)(h 3 (h • X) - z){h • X) 2 
v^X^v^^f(l - (n ■ X) 2 )[l - ^(n-x)-,)^ _ {h ; ^ )2] 



(A.H) 
(A.12) 

(A.13) 
(A.14) 

(A.15) 
(A.16) 
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